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ABSTRACT 


This thesis describes the analysis, method and computer 
programs used to obtain the Fast Fourier Transformation (FFT) 
of an electroencephalogram (EEG) using a small laboratory 
computer like the, PDP Lab 8/E. The EEG power spectrum was 
then computed from this transformation. The information 
contained in this thesis is intended to enable the user to 
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I. INTRODUCTION 


Aeecemeba@al cortex cells discharge electrically, a small 
potential is produced; when this potential is averaged with 
the aeher discharging cells it can be sampled with electrodes 
applied to the scalp. This weak voltage is amplified 100,000 
times by a low noise, high gain amplifier. For simplicity 
and economy, a paper and pen writing oscillograph is used to 
record the changing potential. This record of cortical po- 
tentials is called an electroencephalograph (EEG). 

These waves are very complicated and analyzing them has 
beecneaspreblem singe their discovery. Due to the wave-like 
nature, EEG frequency or spectrum analysis has been the 
subject of much research over the years. Large, high speed 
computers are being used to aid this work today. The pur- 
pose of this study was to develop fast Fourier techniques 
for spectrum analysis of EEGs using a small laboratory com- 
puter lake the PDP Lab 8/E computer system. With the com- 
puter located in the laboratory, near subject and equipment, 
it was hoped that some kind of real time frequency analysis 


could be accomplished. 





II. BACKGROUND 


A. PHYSIOLOGY 

The small potential measured in the EEG originates in 
the outer .1 inch of the brain called the grey matter. The 
working cells of the grey matter are the neurons, numbering 
perhaps to ten billion in all. Very little is known about 
how these cells operate to provide the many functions that 
the brain accomplishes. One well known and supported theory 
by John C. Eccles explains how the EEG waveforms are pro- 
dueed [ReftamelQs)). Bach neuron is assumed to be an independent 
functioning unit, receiving stimuli from outside sensory 
devices and other neurons through dendrite fibers. When a 
particular number of stimuli are received, the neuron fires 
dmpleisesot eclectricity to other neurons via axon fibers. 
Peeetving reverse polarity Signals will cancel the efiect of 
incoming stimuli and prevent the neuron from firing. In 
this way each neuron may be connected to as many as 100 
other neurons. 

Usually many thousands of neurons are stimulated simul- 
taneously. As they reach their excitation points and fire, 
other thousands of neurons receive these new stimuli. Thus 
the collection of stimuli moves through the brain step by 
Seep. wach Stage in the path may be only a thousandth of 
an inch or it may be three or four inches in another part 
Oeeche brain. After firing, each ncuron has a deactivated 


period Of Several seconds, in which it cannot be stimulated 
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to respond. This delay, or refractory period as it is called, 
produces a wave-like quality to the stimuli collection as it 
moves along. Eccles believes that this explains the wave- 
like nature of cortical potentials as they filter out through 
the skull and scalp and are recorded on the EEG. 

Another important aspect of these neural pathways is that 
Eney frequently double back on each other. This causes a 
circular path to form, which provides a self-exciting fea- 
ture to the brain. These closed loops, with fixed distance 
and refractory periods, may account for the predominance of 
certain frequency waves. The prospect of neurons having 
Becdback logps opens Up a whole new level of complexity to 
the brain. In this way each neuron could be viewed as a 
very flexible amplifier, each one operating differently, 


depending on the kind of feed-back it is receiving. 


B. «6©HISTORY 

A short history of the cortical potential must begin 
with Richard Canton, a British physiologist, who reported 
in 1875, "Feeble currents of varying direction pass through 
thesmultiplier when electrodes are placed on two points of 
Pieomexternal surtace (of the brain), or one electrode on the 
grey matter, and one on the surface of the skull." At the 
time, Canton was examining cat, monkey and rabbit brains 
With a sensitive galvanometer using optical magnification. 
It is remarkable that he could distinguish such small po- 


tentials Without the use of electronic amplification, which 





wasn't developed for another fifty years. In 1876 the 
Russian, Danilevskey, discovered a change in the average 
powentval of ithe comtex a.m response to, acoustic stimuli. 
Only recently has this interesting aspect of brain activity 
been reinvestigated. In 1914 Cybulski recorded an apileptic 
seizure in a dog. The early attempts at recording cortical 
potential used photography, however, this was expensive and 
tedious. Pen recorders were not available for this use 
until 1940 [Ref. 8]. 

Recording cortical potentials in man lagged behind their 
demonstration in animals for many years. In man, both elec- 
trodes were attached to the scalp, causing the potentials 
to be greatly attenuated through the skull and scalp. The 
first recordings from the human brain were made in 1924 by 
Hans Berger in Germany. He worked and published much, ex- 
ploring many aspects of the use of the,electroencephalogram. 
Some of these included physiological, neurological, psychia- 
tric and psychological applications. He confirmed Canton's 
findings in animals and established that these waveforms 
change with age, with sensory stimuli and with chemical 
changes in the body. He was the first to record a major 
epileptic seizure in man. Due to the earlier work by 
Kaufman and Cybulski, it became known that epilepsy in 
animals caused abnormal waveforms. 

After Berger's first publication in 1929, many labora- 
Portes Look Up this. work and knowledge in the field grew 


Piel. etitter substantzating Berger's claims, the various 





brain centers began to specialize in different research 
areas: Gibbs and Lennox in epilepsy, the Davises in normal 
subjects, Adrian, Hoagland and Jasper in the physiological 
mechanisms and Walter in brain lesions. The growth of 
electronencephalography has really pavarmeelea = the waeve lop 
ment of the machinery for recording the waves. In the 1930's 
the electronic amplifier gave a tremendous boost to the 
field. In 1948 the transistor allowed the development of 
highly portable equipment. Modern instruments allow 16 
channels to be monitored simultaneously, each of which has 
linear response and accurately reproduces frequencies from 
ecm Weitz emore complete history of electroencephalo- 
GSraphy may be found in Braizier [Ref. 4]. 

One of the many things that Hans Berger discovered about 
the human electroencephalograph was that the most normal 
pattern in adults, with eyes closed, was a frequency between 
6 and 13 Hz whach hewcalled alpha rhythm. These and other 
rhythms seem to originate in the cortex, however, they cease 
if the cortex is surgically disconnected from the thalamus. 
Alpha rhythm is most strong when the subject is relaxed 
with eyes closed. Opening the eyes in normal light signifi- 
Samwely reauces alpha rhythm and similarly, opening the eyes 
in total darkness eae Secciporary recuction. There is 
some evidence [Ref. 18] to support a positive correlation 
between high amplitude alpha rhythms and a passive or in- 
provembed persomality structure as revealed by psychoanal- 
ysis. Similarly, low amplitude alpha rhythms were found 
prevalent iamestroverts [Ref. 18]. 
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Delta activity is between 1 and 3.5 Hz and is generally 
associated with sleep. In 1954 Bremer achieved a strong 
delta state by disconnecting all the sensory pathways. This 
is logical since sleep is most easily induced when sensory 
Stimulation is reduced to a minimum. It must be emphasized 
however, that even though no light, sound, etc. is sensed, 
One does not Sleep unless Ne is “'tired," i.e., delta state 
and sleep are not the same. Theta activity, in the 4 to 7 
Hz range, is usually associated with normal children [Ref. 4]. 
The percentage of theta decreases significantly with age in 
normal subjects. Beta activity, with a frequency of 14 to 
30 Hz, is generally associated with increased anxiety levels. 
The incidence of beta is increased in subjects under-going 
alcohol witndrawai. Beta activity can easiiy be blocked in 
the motor center associated with a particular movement [Ref. 
4). 

The four basic waves, alpha, delta, theta and beta, span 
Elie trequcncy Spectrumeet: cortical potentials. mother pat- 
tern of interest is the mu’rhythm, which is characterized 
Dy a soldat pepeswnenesative half-cyele, rounded positive 
at -Cyeclouandayato liljizeeeMu rhythm is reduced not only 
by overt movement of a limb but also in some cases when 
bie subject is contemplating movement [Ref. 4, 8]. 

Except for beta waves, all the waveforms previously 
@eseribed tend to be reduced by some kind of stimulation. 

In contrast, lambda waves, a saw-toothed waveform, are best 


seenewnen the subject is actively looking at something of 
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interest to him. This waveform usually occurs when the 
subject is scanning a complex pattern or picture. Paradox- 
ically, random waves of the same polarity and shape occur 
in some subjects during light sleep. A thorough analysis 
of the various waveforms in psychologically normal and 
abnormal subjects may be seen in Thorp and Heyman [Ref. 24] 
and Kiloh and Osselton [Ref. 21]. 

Interpretation of the EEG by visual inspection of the 
record depends in a large part of recognizing common pat- 
terns and their association with psychologically normal or 
pathological conditions. Unfortunately these visual methods 
are unable to differentiate between clinical conditions 
having similar EEG's. Several different methods for anal- 
yzing the EEG are being tested. One of the first attempts 
was to get a quantitative grasp on the amplitude of the 
waves. Walter [Ref. 4] joined the peaks and troughs with 
two lines and then measured the average distance between 
Micmemics. tami hemddvent Of amplifiers, electronic inte- 
gration solved this problem. 

In looking for better ways to specify cortical poten- 
tials, many forms of frequency analysis have been tried. 
Boweatly as 1952, Dietsch suggested that these potentials 
might be composed of sine waves and he outlined the procedure 
whereby Fourier analysis could be used to determine the fre- 
quency components present [Ref. 9]. 

ime. tirst frequency analysis of an EEG was accomplished 
Temeecmoinecinnati Coilege of Medicine in 1944, where the fre- 


quencies were determined by manually counting the number of 


LZ 
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complete waves in each second [Ref. 12]. In the late for- 
ties and fifties, much work was stimulated by Gibbs and 
Grass [Ref. 13] when wide-band electronic filters were used 
to separate the different frequencies [Ref. 19, 22]. Filter 
devices are commonly used in this way, even today [Ref. 25]. 
An excellent review of all work on EEG analysis techniques 
up through the fifties may be found in Reference 5S. 

Spectrum analysis of physiological data using the dis- 
crete Fourier transform was accomplished by Halberg in 1961 
(Ref. 14, 23]. Walter was first to use the discrete trans- 
form on cortical potentials in 1963 [Ref. 26, 28]. Most of 
the quantitative frequency work being done on EEG's today 
is with ite aid of large computers in which the discrete 
PoumicTecranstorm provides the conversicn [Ref. 11, 15, 17]. 
vewWe S. Navy is also entering this research area, at the 
Navy Medical Neuropsychiatric Research Unit in San Diego, 
California [Ref. 15]. In 1973 the Fast Fourier Transforn, 

a faster version of the discrete Fourier transform, came 


into use on the EEG [Ref. 2]. 
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IIIT. FOURIER ANALYSIS 


A. FOURIER TRANSFORM 
1. Continuous 
An introduction to Fourier analysis must begin with 
the general form of the Fourier series. It is 
f(t) = Ag + af 


l An cosnut + nZ1 BO Sinnwt (1) 


where w = 27/T and T is the period. Any waveform can be 
represented by determining the coefficients Ao» AL and BL 
and summing all the terms of this series. The coefficients 
give the relative contribution of each frequency that is 
present in the signal. These are harmonics or integer 
multiples of the fundamental frequency 1/T. Fourier de- 
veloped a set of equations from (1) and solved them simul- 


Gameously. ihe result 


> 
O 
N 
‘oe | Lee 
Hy 
on 
ct 
Nena 
Ou 
ct 


(2) 
O 
5 T 
A= f f(t) cosnwtdt (3) 
2 1a | 
BO =F £ fGeies i nnwedt (4) 
O 


provides a method for finding the harmonics when the function 
is known. 
Zee Discrete 
Frequently, however, the complete function is not 


known and only a discrete number N of samples over a finite 
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period T is known. In this situation, an approximation to 


the continuous function may be written 


n n 
tte A + 2H] A, cosmwt + 202 BL Sinmwt (5) 


O 1 


where n = N/2 for N even and n = N+1/2 for N odd. The ex- 


pressions for the coefficients now become 


pede ae) 2mmt 
An ~ WV re-n f(t) COS ( VY ) (6) 
a f(t) sin ( 2nmt ) (4 
m N r=-n N : 


These functions are easily computed from the set of 
Sample data points; however, the time required is sometimes 
prohibitive. It will be noted that this procedure requires 
N? trigonometric operations, since each of N coefficients is 
computed from N terms of a series. These computations, in 
addition to numerous other arithmetic operations, require 
considerable computer time, especially on small laboratory 


computers like the PDP Lab 8/E. 


B. PROBLEM AREAS 

There are three problems commonly encountered when using 
the discrete Fourier transform. These are: aliasing, leak- 
age and the picket-fence effect. Aliasing occurs when a 
heer s than expected frequency is sampled, but appears in 
the spectrum as a lower frequency. Suppose an unknown sig- 
nal is being sampled five times at a rate of 100 Hz. If 


an 80 Hz component is present, it will appear in the spectrum 


Ie 





as 20 Hz and there is no way of distinguishing this pseudo- 
frequency from the real one. Nyquist proposed the solution, 
in that, the sampling rate should be twice the highest fre- 
quency present [Ref. 16]. This can’be accomplished in an 
unknown signal by passing it through a lowpass filter set 

at half the sampling rate. A lowpass filter was used to 
prevent aliasing in this study. 

The problem of leakage occurs in all Fourier analysis 
using finite sampling records. The record, or the total 
length of the sample, can be viewed as a window. When this 
window is transformed, without the sample data present, the 
resulting spectrum contains not only the fundamental fre- 
quency but also spurious peaks called sidelobes. When the 
window with its sampie is transformed, these sidelobes tend 
Eomredk into the spectrum. Because of the relatively small 
Sidelobes, in the frequency range of interest, leakage was 
not a problem in this study. 

The packee snes Cimeme 1S a problem when a pure sinue 
soid occurs between two adjacent harmonics. If it is ex- 
actly halfway between computed harmonics, the signal is 
reduced to .637 of the original in both harmonics. The 
worst situation occurs when there is a frequency between 
each pair of adjacent harmonics. In this case an erroneous 
Tippie, or picket-fence, appears in the spectrum. Pure 
Sinusoids were assumed not to be present in EEG's and thus 
the picket-fence effect was ignored. Bergland discusses 
these problems So as the mathematics associated with 
them in Reference 3. 
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IV. FAST FOURIER ANALYSIS 


A. FAST FOURIER TRANSFORM 

In 1965 James Cooley and John Tukey published an impor- 
tant paper that described a much faster method to compute 
the discrete Fourier transformation, called the Fast Fourier 
Transform (FFT) [Ref. 7]. The procedure adopted for use on 
high speed computers has become known as the Cooley-Tukey 
Algorithm. A general outline of the proof of this algorithm 
follows. 

A time series Xe containing N even samples, can be par- 
titioned into two functions i and Zi. each containing N/2 


peaints. Let io Demeomposecdmot the ecVvell-numpered Ponts 


(X_sXo9-+- Xu) and Z.. the odd-numbered points (Xq5Xz---Xy_y2)- 


Equations (6) and (7) may be combined using complex number 


notation, so that 


n-1 
2 mr 
am 720 xX. W (8) 
where W = ge and j3 = v-1l. Since . and Z. Aa) Bl SNe oye 


N/2 points, each has a discrete Fourier transform whose 


coefficients are 


M2 1 ae 
Bee. 2 om nal (9) 
N/2-1 
_ m(2r+1) 
with m = OREM). ws y7 1. Combaninegsedquations (9) and (10), to 


get odd and even points, gives 


lif 





N/2-1 
_ mr m(2rt+1) 
Peele iY Wo + 2 ] 
N/2- N/2-1 
= mr m mr 
ore tT WE 2 WN 
A = B + Wc (11) 
A +N/2=B - w" Cae (12) 


The details of this proof will be found in Reference 6. 

pe x. is an even power of two, then ie and Z. also have 
an even power of two points and each can be partitioned again 
and again until only one point remains in each set. The 
Fourier transform of a single point is just that point again; 
the entire transformation can be accomplished using equations 
(11) and (12) repeatedly at each partition. The requirement 
that the number of points N be a power of two assures that 
ee ce Tr log, N partitions there will be only one point in each 
Bet In this way, the entire set of points can be translated 
without carrying out the manipulations of equations (6) and 
(7). Since the FFT is derived from the discrete Fourier 


transform, the FFT maintains the same precision as before. 


B. FLOWPATH 

An example of the FFT for N = 8 will be illustrative. 
Figure 1 shows the flow diagram and the configuration after 
the first partition. Using Tukey's notation, each dot rep- 
Pee @esed Variable and the factors entering the dot are 


additive. When a multiplication is required before adding, 
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twesdUantity will be beside the factor to be multiplied. 
The set of A. for j = 0,...,7 is computed from equations 
weeyeond (12). The quantity near the bottom right dot is 
A. = B, + w’C.. Figure 2 shows all quantities after two 
partitions and Figure 3 shows’ the complete transformation 
broken down into only multiplication and addition operations. 
It may be noted in Figure 3 that the data points must be 
rearranged according to the bit reversal rule before the 
algorithm can be carried out. Bit reversal is accomplished 
by representing the data point position number in binary 
and reversing the bits about the center. Position (3) is 
Peeeeecim Owl jeand 1s reversed to (1 1 0), which is (6). 
Thus (0) stays a (0), (1) becomes a (4) and so on. More 


information about FFT may be found in References 3 and 6. 


C. PROGRAM 
The computer program FAST was written to accomplish the 
Heleon them’ DP Lab 8/E computer system. A complete listing 
of FAST may be found in Appendix A. The program is written 
in BASIC, one of the commonly used languages on the PDP Lab 
8/E. FAST receives data points as inputs and outputs the 
weurier coefficients. 
timoneeustOuoperate the FFT, read FAST into core using 
normal methods, and Benepe RUN. The program will type 
ENTER DATA POINTS 
U 
Picmeomilto a Wait loop. Type in 64 data points in their 


proper. order. The number of data points may be changed to 
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any power of two, but cannot exceed 64 due to core limita- 
tions, which is 8,000 words in the PDP Lab 8/E. FAST will 
translate the data points and type the Fourier coefficients. 
The time required to convert 64 points is approximately 25 
seconds. The power associated with each harmonic can be 
obtained by squaring and adding the coefficients at each 


frequency. 


D. EXAMPLE 

As an example to illustrate the use of FAST, the square 
wave shown in Figure 4 was translated and the results com- 
pared with the coefficients from the continuous Fourier 
transform: Table I shows the result of a FAST translation 
Oit,ehne square wave along with Che coerficients. To deter- 
mine the AL coefficients using the continuous transform, 
equations (2) and (3) are integrated from -T/2 to T/2. 
SIcce mi t)) 1S dneodd function, thie aan terms are ail zero: 
Using equation (4), however, and integrating over the same 
range, it is easily shown that 


n nt 


To compare with FAST, the same 32 harmonics were computed 
using equation (13)(see Table II). The answers from FASI 
are within two decimal places of the answers computed from 


the continuous Fourier transform. 
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Figure S. EEG Wave and Timing Pulse. 
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HARMONIC A(N) B(N) 


0 0 0 

1 0 1A Ata ay 

2 0 0 

3 0 ~4213407 
4 0° 0 

5 (0 - 2495139 
6 0 0 

7 0 ~ 1746758 
8 0 0 

2 0 ~1321451 
10 0 0 

11 0 ~ 1042749 
12 0 0 

13 0 .08427145 
14 0 0 

15 get .06895806 
16 0 0 

17 0 .05664668 
18 0 0 

19 0 -04635313 
20 0 0 

ZA 0 -03746103 
ae 0 0 

23 0 .02956024 
24 0 0 

Z5 0 .02236283 
26 0 0 

if 0 .01565539 
28 0 0 

Zo 0 Sey Cagis yop 10S) 
50 0 0 

31 0 SUT UZOSE= 0S 
READY. 


Rabplewls FAST Coefficients Translated from the Square Wave 
of Figure 4. 
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RUN 


HARMONIC B(N) 
0 5.50 FOO bao 7 
il 275241 
2 “A. 74 SLO TE-09 
5 ~4244135 
4 a 2 Suelo LU 
> -2546481 
6 = Tyo o Vote 09 
7 wloleogls 
8 Sky high JAS HoH otters OS 
9 ~1414712 
10 -9.486382E-10 
ili whls7 491 
£2 0S 5155-10 
3 .09794159 
14 DOr 967 E- ko 
15 .08488271 
16 | Boromoo coe 10 
Ly .07489651 
18 one U2 Ze 10 
19 SUG70 1207 
20 oA ao B= 10 
21 .06063051 
22 -4,311992E-10 
ZS BUS 905029 
24 =e 5Z2O59E- 10 
As -05092963 
26 -3.648608E-10 
a .04715706 
28 aye Sis S/S Nein 
29 .04390485 
30 ESO Zi 7 E10 
Sl .04107228 
READY 

Table II. Continuous Fourier Transform of BO Ceerticients 


fermesduare Wave of Figure 4. 
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V. SPECTRUM ANALYSIS 


he oOFECTRUM 

In order to do frequency analysis of EEG waveforms, a 
computer program called SPECTRUM was written which incor- 
porated FAST (Appendix A). The problems encountered were 
many. First, a relatively noise-free signal had to be pre- 
sented to the PDP Lab 8/E analog-to-digital (A-D) converter. 
Then the signal had to be sampled at the right rate and 
duration. Finally, the data had to be delivered to FAST in 
a compatible format for translation. The power spectrum was 
chosen as the output format since it has more intuitive ap- 
peal than some other measure such as voltage or current. 
According to Parseval's theorem, the power at each frequency 
can be computed by squaring the two coefficients and adding 
them together. This was done and the power spectrum is dis- 
mrayed on tite teletype as percent of total power at each 


harmonic. 


B. APPARATUS 

The Beckman type RM Dynograph Recorder was used to pick- 
moeand amplify the cortical potentials. Using sodium chlo- 
mueeeDaste, Silver electrodes were applied to the subject's 
scalp at positions C, and Tp. A standard Beckman A-C 
Seupler was used with the TIME CONSTANT switch set to .05 
and HIGH FREQ switch on 3. Depending on the signal strength, 


the preamplifier was set to either 1 mv/cm or .5 mv/cm and 
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tne power amplifier set for a multiplication of .02. This 
adjustment was made in order to maintain a relatively con- 
Beart Signal level from the power amplifier. The FILTER 
switch was set on 3, which is a 60 Hz lowpass filter. 

The EEG was recorded on one track of a Hewlett-Packard, 
four track, instrumentation recorder on 1.0 mil magnetic 
tape at 3-3/4 in/s tape speed. A seven ms timing pulse was 
recorded on another track from an Ortec pulSe generator. 
The timing pulse made it possible for the computer to find 
a particular place on the tape. The two tracks were then 
Eomnected toutne A-D converter of the PDP Lab 8/E computer. 
As an alternative hookup, the A-D converter can be monneetca 
directly to the Beckman and sample the EEG as the subject 


proawees it. 


C. SAMPLING RATE 

PEC imeunicmDMonpersSional on channel, 1 jof the A-D converter, 
SPECTRUM samples and translates it to the power spectrum. 
According to Nyquist's rule, the sampling rate must be twice 
the maximum frequency. The maximum frequency passed by the 
Beckman filter is 60 Hz, so the sampling rate must be at 
PaierUed?; weeiieorder to give a little extra margin, 128 Hz 
was chosen so that the harmonics would come out in even 
multiples of two. FAST can only convert 64 points and sam- 
pling at 128 Hz, the period (T) of the sample is % second. 
Phesteciprocal of the period (1/T) determines the increment 
of the harmonic. In this case 1/T = 2, thus each harmonic 


Pmexaetiy 2 Hz apart. 
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D. EXAMPLE 

After loading SPECTRUM and typing RUN, the program sam- 
ples, converts the data, and types the power spectrum. 
Table III shows an example of one SPECTRUM run, taken from 
a subject who was relaxing with his eyes closed. The pro- 
gram starts sampling when it receives a numbered timing 
pulse on channel three. The EEG waveform from which this 
example was taken, along with its timing pulse, is shown 
in Figure 5. SPECTRUM uses the timing pulses to locate the 


desired sample on the tape. 


Be. EXPERIMENT 

By choosing the proper segment with the timing marks and 
replaying the tape 20 times, SPECTRUM was used to scan a 
.length of tape arbitrarily set for 10 seconds and produce 
the power spectrum for each half second record. The results 
of this experiment were reduced to the graph of Table IV. 
Meeplots frequency against time, with each number represent- 
Migeune Mercentage Of total power for a given frequency and 
time. If the table is viewed as a three dimensional graph, 
where percentage power is shown coming out of the page, it 
1s easy to visualize a ridge running up the 10 Hz line. The 
subject mae again resting with eyes closed, so the prepon- 
derance of alpha frequency (8, 10, 12 Hz) was understandable. 
The frequency also appeared to come in bursts. It was first 
on for about three seconds, then reduced for one second, then 


back on for a second, and so on for the length of the tape. 
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RUN 


FREQUENCY . POWER 
HZ % 
0 Ge 2 Los 
1 1915/7298 
Us .07388068 
) 4.684211 
4 5.177504 
5 a7 29045 
6 49.26551 
SLO 7 
8 Heh lhe ltes Ie 
a Sa 954.7 5 
10 . 7096686 
11 O07 305 
eZ res) See) Us 
13 70 O01 
14 -6030759 
15 soo 0055 
16 eo Lo 
i 22905483 
Se ~ 3325844 
19 i299 O47 5 
20 SDT OS MOL 
21 Ey, 1572 (ors) 8 
Age BAS Sees 
25 ae S10 9 
24 ~1704812 
Pigs: .1763688 
26 ~1824439 
Ze 7 los Legs 
28 eel Honelge: syal 
VS) ~1526896 
30 . 1463384 
31 ~ 148733 
READY . 


acm eee aA lypical SPECTRUM Rum from a Half Second Sample 
Sed BEG . 
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Table IV. 


A 6 CuO eee 14 toe 1S. 20 22 24 26 28 350 


FREQUENCY (Hz) 


Power Spectrum derived from an EEG Ower a Penaod of 
Loimsceonas, veae1 Nimmer Represents the Percentage 
Power aisecach Erequemey for the Half Second Sampled. 
ZeTros are ,omitted. 
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Sometimes this alpha burst effect can be detected right on 
faite ooG. 

There were two places where strong delta frequency was 
meesent, Prepresented on the graph by the 2 Hz line, delta 
was present in the beginning and again around the five 
second point. This normally indicates a very light stage 
of sleep, where the subject is drifting in and out of con- 
sciousness. Under questioning the subject indicated that 
he could have been asleep several times during the recording 
session. The graph also shows the rapidity with which fre- 
memey chamocominert. 2, 11, 15 and 27]. One example of this 
is from the first half second to the one second point, where 
the 10 Hz power changed from 3% to 39%. Rather than changing 
smoothly, it was as if one frequency stopped and another 


Started. 


me)6UNOISE 

Either because of a noisy amplifier or possibly a dirty 
tape, it may become necessary at some time to analyze a 
noisy signal. This will become known when repeated spectra 
from the same point show a variance. In order to reduce 
this variance, a method of averaging was devised. Instead 
of Sau) bee Saliceei2, a mew rate of I280 Hz was used and 
the sample means were computed every 10 points. The sample 
period was kept at % second, so a total of 64 sample means 
were computed each run. This reduced the variance by a fac- 
tor of VI0 and providing the signal’ was not too noisy, the 


results were satisfactory. Thus to obtain more information 
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from the noisy EEG, it may become necessary to sample more 
frequently than according to Nyquist's law. 

Increasing the sampling rate is not without its short- 
Seemines. Firstly, it takes longer to process the added data. 
Changing to a 10 point average increases the computing time 
om 25 to 40 seconds. Secondly, part of the data must be 
stored while the averages are being computed. The new pro- 
gram requires about 200 more storage positions, which is at 
the maximum capability of the PDP Lab 8/E with 8,000 word 
memory. Most spectral analysis work being done on EEG's 


today use a 200 Hz sampling rate [Ref. 2, 15]. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A real time analyzer for the EEG was one of the goals 
of this study. This has only partly been realized since 
SPECTRUM requires 25 seconds to convert a half second sam- 
micas oimnce the EEG frequency pattern changes SOmigcod a |yaeelae 
1s necessary to convert more than a half second; perhaps as 
Ruch as 10 seconds in order to get a good representation of 
the changing nature of the wave. This means. that a much 
larger computer is needed to do real time spectral analysis 
ot EBEG'S. 

As was shown using off-line techniques, SPECTRUM is a 
useful tool for spectral work. Using an instrumentation 
recorder in conjunction with some kind of pulsing device 
moretiming, many different forms of frequency work may be 
accomplished such as autocorrelation, cross spectral analysis 
and digital filtering. SPECTRUM can also be used to compute 
power spectra from other kinds of continuous waveforms such 
as electrocardiograms and electromyograms. 

iiewavutnor believes that this area of study is an impor- 
tant one and should be extended in the future; however, 
there are iiiecterons, (he major limitation to this system 
meetaes in the lack of adequate data storage, both in the 
small core size, which is about 4,000 words with BASIC in 
core, and in the lack of a random access device such as a 
disk pack. Of these two, the disk pack would be most useful 


for spectrum analysis work since large amounts of data from 
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the A-D converter could be stored rapidly and processed as 
required. The flexibility of the entire system would be 

considerably enhanced because the disk pack could also be 
used to break up large programs which are now too big for 


the computer and then process them in a segmented fashion. 
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APPENDIX A 


SPECTRUM 
1 N=64 
Z N1=63 


fee DIM F(1,63) ,W(1,31),T(1) 
8 SET RATE 6,558 
10 DIM P1(300) 
15 REAL TIME P1,1,0,900 
18 IF ADC(3)<=.9GO TO 18 
20 ACCEPT 
30 FOR I=0 TO 63 
35 FOR J=1 TO 14 
40 F(0,1)=F(0,1)+ADB(0) 
J 


45 NEXT 
50 NEXT I 
oo PRINT 


100 FOR I=0 TO 31 

110 READ F(1,1) 

115 F(1,63-1I)=63-F(1,1) 
120 NEXT I 

158 FOR I=0 TO Nl 

162 IF F(i,i)=i THEN 169 
163 X=F(0,I) 

164 J=F(1,I) 

165 F(0,I)=F(0,J) 

166 F(0,J)=X 

167 F(1,1)=F(1,J) 

168 F(1,J)=J 

169 NEXT I 

194 FOR I=0 TO Nl 

195 F(1,I)=0 

196 NEXT I 

210 W(0,0)=1 

220 W(0,1)=COS(6.283185/N) 
230 W(1,1)=-SIN(6.283185/N) 
240 IF N<=4 THEN 295 

250 FOR I=2 TO N/2-1 

260 I1=I-1 

270 W(0,I)=W(0,1)*W(0.,I11)-W(1,1)*W(1,I1) 
280 W(1,I)=W(1,1)*W(0,11)+W(0,1)*W(1,11) 
290 NEXT I 

310 J=1 

320 FOR M=1 TO N 

330 J=J*2 

340 IF J>=N THEN 370 

350 NEXT M 

370 FOR K=1 TO M 

380 J2=2t(K-1) 

390 I2=N/(J2*2) 
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400 FOR I=0 TO 12-1 

401 T(0)=F(0,Y) 

410 FOR J=0 TO J2-1 

500 X=N*I/12+J 

510 Y=X+J2 

520 Z=J*12 

530 T(0)=W(0,Z)*F(0,Y)-W(1,Z)*F(1,Y) 
540 T(1)=W(1,Z)*F(0,Y)+W(0,Z)*F(1,Y) 
550 F(0,Y)=F(0,X)-T(0) 

Geo F(i,Y)=F(1,X)-T(1) 

570 F(0,X)=F(0,X)+T(0) 

Boo F@l,X)=F(1,X)+T(1) 

590 NEXT J 

600 NEXT I 

610 NEXT K 

690 A=0 

700 FOR I=0 TO N/2-1 

710 F(0,1)=(F(0,1)+2)+(F(1,1)+2) 

712 A= A*F(O, I) 

720 NEXT 

728 bee " FREQUENCY POWER" 
730 FOR I=0 TO N/2-1 

732 A2=F(0,I)/A 

750 PRINT 1*2,A2 

760 NEXT I 

800 FOR I=0 TO Ni 

810 F(0,I)=0 

820 NEXT I 

900 DATA 0,32,16,48,8,40,24,56,4,36,20,52,12,44,28,60,2 
901 DATA 34,18,50,10,42,26,58,6,38,22,54,14,46,30,62,1 
998 RESTORE 

999 GO TO 1 


al N=64 

g N1=63 

5 ites, 65) ,N(1,31),T C1) 
30 FOR I=0 TO Ni 

mo INPUT F(0,I) 


49 F(1,I)= 
50 NEXT I 
58 PRINT 


100 FOR I=0 TO 63 

io READ F(1,1) 

120 NEXT I 

130 FOR I=0 TO N1 

140 I1=1*64/N 

150 F(1,1I)=F(1,1I1) 

155 NEXT I 

158 FOR I=0 TO N1 

162 IF F(1,1)=I THEN 169 


36 





163 X=F(0,I) 

164 J=F(1,I1) 

165 F(0,I)=F(0,J) 

166 F(0,J)=X 

167 F(1,I)=F(1,J) 

168 F(1,J)=J 

169 NEXT I 

194 FOR I=0 TO Nl 

195 F(1,I)=0 

196 NEXT I 

210 W(0,0)=1 

220 W(0,1)=COS(6.283185/N) 

230 W(1,I)=-SIN(6.283185/N) 

240 IF N<=4 THEN 295 

og FOR 1=2 TO)N/2-1 

260 I1=I-1 

270 W(0,I)=W(0,1)*W(0,11)-W(1,1)*W(1,I1) 
280 W(1,I)=W(1,1)*W(0,11)+W(0,1)*W(1, 11) 
290 NEXT I 


510 J=l1 

320 FOR M=1 TON 

S50 J=J*Z 

340 IF J>=N THEN 370 
350 NEXT M 


370 FOR K=1 TOM 

380 J2=2+(K-1) 

390 I2=N/(J2*2) 

400 FOR I=0 TO I2-1 

401 T(0)=F(0,Y) 

410 FOR J=0 TO J2-1 

500 X=N*I/I2+J 

510 Y=X+J2 

520 Z=J*I2 

530 T(O)=W(0,Z)*F(0,Y)-W(1,Z)*F(1,Y) 
S540 T(1)=W(1,Z)*F(0,Y)+W(0,Z)*F(1,Y) 
550 F(0,Y)=F(0,X)-T(0) 

560 F(1,Y)=F(1,X)-T(1) 

570 F(0,X)=F(0,X)+T(0) 

580 F(1,X)=F(1,X)+T(1) 

590 NEXT J 

600 NEXT I 

610 NEXT K 

698 A=0 

699 PRINT "HARMONIC A(N) B(N)" 
700 FOR I=0 TO N/2-1 

mmo PRINT I,F(0,1)/32,F(1,1}/32 

720 NEXT I 

730 FOR I=0 TO N/2-1 

736 V1=2*I/N 

738 V2=F(0,I)/A 

739 IF V2<1.000000E-03GO TO 760 

740 PRINT I*2:V2 

760 NEXT I 

761 STOP 
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800 FOR I=0 TO N1 

810 F(0,1I)=0 

820 NEXT I 

900 DATA 0,32,16,48,8,40,24,56,4,36,20,52,12,44,28,60,2 
901 DATA 34,18,50,10,42,26,58,6,38,22,54, 14,46, 30,62,1,33 
Boz DATA 17,49,9,41,25,57,5,37,21,55,13,45,29,61,3, 35,19 
5 DATA 51,11,45.27,59,7,39,23,55,15,47,31,63 

998 RESTORE 

999 GO TO 1 


Sore) 





EO. 


aeall 


a2. 


BIBLIOGRAPHY 


Anand, B. K., Chhina, G. S. and Singh, B., "Some Aspects 
of Electroencephalographic Studies in Yogis," 


ble GunecieciemeimmemNcuUropmys10l.,, W961, 13: 452-456. 


Banquet, J. P., "Spectral Analysis of the EEG in Medi- 


tatiom Electroeneeph. clin. Neurophysiol., 1973, 
5 5 eaeeses = 15 1h, 


Bergland, G. D., "A Guided Tour of the Fast Fourier 
Creams Toriy weber opeetrum yy 1969, 6" 41-52. 


Brazvet we Wucuulectrical Activity of the Brainy 
Williams and Wilkins Company, 1968. 


Burch, N. R., "Automatic Analysis of the EEG: A Review 
Heke sitmeatlOne Oot Systems," Electroenceph. clin. 
Neunopiy si ole 59, 11 semei27-854. 


Cemeran wien. , Cooley, J. Ws, Favin, D. L., Helms, 
feo, Kacneloerk. A., Land, W. W., Malling, G. C., 
Nelson, DS E,,@Rader, C. M. and Welch, P. D., 
minat 1s the Fast Fourier Transform?,'* IEEE Trans. 
Audio and Electroacos., 1967, AU-15: 45-55. 


Cooley, J..W. and Tukeys J. W., #An Algorithm for the 


Machine Calculation of Complex Fourier Series," 


Natvematlcamorecomputation, 1965, 19:0 297-301. 


Cooper wee, Osselton, J. W. and Shaw, J. C., EEG Tech- 
MolcGeyemoputterworth and Co., Ltd.,)1969- 


Dietsch, G., "Fourier-Analyse von Elektrencephalogrammen 
Giememienm=schen,'' bPilugers Archives, 1952, 230: L06-112. 


BCGlCrn) mem he Physiology of Gmagination,” Scien- 
tific American,-1958. 


Elazar, Z. and Adey, W. R., “Spectral Analysis of Low 
Frequency Components in the Electrical Activity of 
the Hippocampus During Learning," Electroenceph. 


Smmpmeencunopnysiol., 1967, 257°" 225-240. 


ree ome. L., Romano, J., Ferris,eE. B., Webb, J. P. 
and Stevens, C. D., "A Simple Method of Determining 
Frequency Spectrums in the Electroencephalogram," 


mmepives Of Neurology and Psychiatry, 1944, 51: 
lead -146. | 


a9 





MS. 


14. 


VS. 


16. 


7 


ES: 


Lee 


BAO. 


21. 


Oe. 


o5. 


24. 


ES . 


Gibbs, F. A., and Grass, A. M., ''Frequency Analysis of 
Powe “cGhenececomaga,, January 253, pp. 152. 


Halberg, F. and Panofsky, H., "Thermo-Variance Spectra: 
Method and Clinical Illustrations," Experimental 


MedicalmoumcenyyeelSOl, 19: Z264=3509. 


~Hord,, D.. Je iacene fe. and Johnson, Ls C.... "Intensity 


and Coherence Contours During Self-Regulated High 
Alpha Activity," Report Number 17-54, Navy Medical 
Neurophychiatric Research Unit, 1971: 


Jenkims More tieedncandtts. DD. G. > Spectral Analysis and 
lismopolucatone, molden-Day, Incs; 1968. 


JONSON Ubi.) Nattch, P., Nute, C. and Austin, 
M., "Spectral Analysis of the EEG of Dominant and 
Non-Dominat Alpha Subjects During Waking and Sleep- 


i eee eonocnecern. clin, Neurophysiol., 1969, 26: 
Sig lkrow, On 


NKCChe Oma Ned, Eo. Applied Physiology, Oxtord 
Undvemsary Press, 1961, 


KeontardametienA., Rabinovitch, M. S. and Fister, W. P., 
"The Use of Frequency Analysis in the Interpretation 
of the EEG of Patients with Psychological Disorders," 


blic@mmeenecenn. Clin. Neurophysiol: ; 1955, 7: 29-38. 


Kennard, M. A. and Schwartzman, A. E., "A Longitudinal 
Study of EEG Frequency Patterns in Mental Hospital 
Fattcmesmana Normal Controls,” Electroenceph. clin. 
Neurem@ietol.,, 1957, 10: 263-274. 


Ki lon wameocnedand Osselton, J. W., Clinical Electroen- 
Cepligeaemapny, Butterworth and Co.;, Ltd., 1966. 


Kettichpwreeand Margerison, J. H., “Présentation of 
Data From Low Frequency Analysis to Illustrate 


Scmlcmmenanges in the EEG)” Ellectroenceph. clin. 
NEWmopmystol., 1961, 13: 606=61T- 


Panofsky, H. and Halberg, F., "Thermo-Variance Spectra: 
Simplified Computational Example and Other Methodo- 


tO eee xperimental Medical Surgery, 1961, 19: 
aor OG 


iemoewneG. and Heyman, A., Principals of Electroen- 
eopmalography, Duke University Medical Center Press, 
1966. 


Veldavka, J., Watousek, M. and Roubicek, J., "Mental 
Arithmetic and Eye Opening. An EEG Frequency Anal- 


PeiemaideGoak Study,  Ellectroencepn. Clin. Neuro- 
Piero, 1967, 22: 174-176. 


40 





woeeeWalter, D. O., “Spectral Analysis for EEG: Mathemati- 
cal Determination of Neurophysiological Rationships 


from Records of Limited Duration," Experimental 
Neurology, 1963, 8: 155-181. 


femevatcer, 0. O. and Adey, W. R., "Spectral Analysis of 
PeG@meRecordea During Learning in the Cat, Before 


and After Subthalamic Lesions," Experimental Neur- 
oleae, (2 481-501. 


28. Walter, D. O., Rhodes, J. M., Brown, D. and Adey, W. R 
"Comprehensive Spectral Analysis of Human EEG Gen- 


Pltmomm OsteriOr Cerebral Regions,” Electroenceph. 
Cuiimectnophysioiy 1966, 20: 224-237. 


ous 


41 





INITIAL DISTRIBUTION LIST 


No. 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Department Chairman, Code 55 
Department of Operations Research 
and Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93940 


Assoc Professor G. K. Poock, Code 55 

Department of Operations Research 
and Administrative Sciences 

Naval Postgraduate School 

Nowtete meat ttoria 93940 


Assoc Professor T. A. Wyatt, Code 55 

Department of Operations Research 
and Administrative Sciences 

Naval Postgraduate School 

Monterey, California 93940 


Assoc Fromessor R. Bedow, Code 55 
Department of Operations Research 
and Administrative Sciences 
Naval Postgraduate School 
MonveveyewGalttorniia 93940 


Chief of Naval Personnel 
(Pers 11b) 

Department of the Navy 
masnameceom, D. C. 203570 


LT Lawrence M. Gorham, USN 


507 Carmel Avenue 
Marina, California 93933 


42 


CCpiles 














thesG5764 
EEG frequency analysis on the PDP Lab 8/ 


iii 





